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I use the numerical values of the generalised Glaisher-Kinkerlin- 
Bendersky (GKB) constants to give numerical values for the deriva¬ 
tives of the Hurwitz (^-function at negative integers, rather than the 
other way round. I point out that both Glaisher’s numerical approach 
and Bendersky’s recursion for the generalised gamma function were 
anticipated by Jeffery in 1862 who gave the value of the second con¬ 
stant as an example. I therefore propose that GKB become GKBJ. 
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1. Introduction 


In many calcnlations of the effective action, long combinations of derivatives of 
Hnrwitz (^-fnnctions occnr and often it can be more nsefnl to give their n um erical 
valnes. Althongh this is a fairly standard topic, some farther discnssion might be 
of interest. One conld of conrse use a built-in algorithm in a suitable CAS, but 
if this is not available, a direct analytical calculation is necessary, preferably not 
involving quadrature. Suitable forms have been found by Elizalde, [1], some time 
ago. My calculation will incidentally produce the same (asymptotic) expression but 
by a different route. 

2. The basics 

I approach the formula in a possibly circular way allowing me to introduce 
some other notions and, later, a few historical remarks. I start with the basic 
relation between and ^^(0,1), the Riemann (^-function, with (initially) 

vj an integer, 

C'(-fc,«;)-C'(-fc)-logr,(«;), (1) 

dropping the ‘H’. Here rj[.(r(;) is Bendersky’s generalised gamma function dehned 
(initially) by the sum 

W 

^ogTkiw + 1) = Z logm , w G Z. (2) 

m=l 

There are two attitudes one can take towards (1). These are, in essence, 
whether it is to be read from right to left or left to right. In the hrst case, which 
is the one I mostly adopt here, the right-hand side, dehned as the explicit sum 
(2), is analysed independently using hnite calculus summation techniques leading 
to calculable quantities in which w can be considered continuous. This is the proce¬ 
dure of Bendersky, [2], and Jeffery, [3]. (They do not consider the relation with the 
(^-function.) Having an expression for the generalised gamma function then allows 
information about the ^-function to be extracted from (1). 

On the other hand if the C^function is assumed known then T;. can be taken 
as defined by (1) and its properties determined therefrom. I will not consider this 
point of view much here but it has the advantages of rapidity and elegance. An 
example will appear later. 
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I now invoke, withont worrying too mnch abont where it conies from, Adam¬ 
chik’s formnla for ('{—k), [4], 

C'{-k) = !^f^-logA„ (3) 

where are Bendersky’s generalisation of the Glaisher-Kinkelin constants dehned 
by the asymptotic behavionr. 


logA, = L, = J^logr^(w; + 1) ^_i„dependant ' 


( 4 ) 


is a harmonic nnmber. 

One writes, 

logrfc(M; + l) =Lfc + Afc(M; + l) (5) 

where Ai^{w) is a fnnction with the appropriate limiting behavionr. Bendersky de¬ 
rives an explicit asymptotic (divergent) series from the Enler-Maclanrin snmmation 
formnla applied to logr;.(r(; + 1). It is not necessary to write it ont yet. 

Combining (1), (3) and (5) trivially yields the reqnired compntable formnla. 


C{-k,w) 


k + 1 

^k+i 

k + 1 


- Lk + logVkiw) 

+ Afc(w), 


( 6 ) 


which gives the asymptotic expansion obtained by Elizalde directly from integral 
forms of the ^-fnnction. 

I look npon (3) as a means of Ending ('{—k) from the A^, rather than the 
other way ronnd, as is more nsnal. As a check of the nnmbers, if k is even ('{—k) 
can be transformed into C(2/j + 1) the valnes of which are readily available, to high 
accnracy.^ 

Nnmerically one can treat the two lines in (6) separately. In the top line 
is to be calculated by trial from (4) using the known asymptotic form of A. (The 
method used by Jeffery, Glaisher and Bendersky). Less accurately, one can just 
substitute the asymptotic form directly into the second line. Of course, using an 
asymptotic form is not the best numerical procedure, but an accuracy of more than 
10 places is easily attained this way. (The isolated case of re = 1 is not possible in 
the second approach.) 

^Or say Wolfram Alpha can be employed. 
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The more accurate computation follows from the hrst line in (6) where is 
found from (4) using the asymptotic form for Ai^{w + 1) with a suitable choice for 
w. logrj!.(r(;) is determined ‘exactly’ by (2). This method is therefore restricted to 
w integral. 

If w is not integral, one has to use the second line in (6). If rc is small a direct 
application is not possible but one can translate rc to a large enough value by adding 
an integer and then employing the fundamental property of the generalised Gamma 
functions, [2], valid for any x. 


rfc(a; + l) 


( 7 ) 


several times, if necessary 


3. The details 


I will now hll out the previous discussion by outlining Bendersky’s approach. 
I do this because his paper is not widely recognised. To begin with, it is helpful 
to write down the explicit expression for the asymptotic series, A^i., as given by 
Bendersky. One particular form of this is, and I write it out exactly as in [2] p.276, 
(with a few misprints corrected). 


Afc+i(a: + 1) = 


X 


k+2 


k + 2 

fc-i 


log a: 


X 


k+2 


+ (k + iy.J2 


s, 


r+1 


(/c + 2)2 

xk+\—r 


+ -X 


k+1 


loga:+ 


(r + 1)! (/c + 1 - r)! 


log a; + 


1 


+ ^loga: + 

OO 

+ (k + 


/c + 1 k 

1 n B 


1 

+ ^ + ...+ 


k + 2 — r 


+ 




S = 1 


Bk+i+s (’S~2)! 
(/c + l + s)! x‘^~^ 


( 8 ) 


which is easily coded. 

Bendersky derives this, at general k, from an Euler-Maclaurin summation, 
after some amalgamation of terms. The details are of no immediate concern. 

It is always illuminating to have some specihc examples before one’s eyes and 
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I give the following, 


. , , x{x + l){2x + l) , x^ X 1 1 

KJx + 1) = —- - logo:-^-^- 5 - 

^ 1.2.3 ^ 9 12 360 a: 7560a:3 

1 


K^{x + 1) = 


x{x + 1) 1 \ x^ 1 

L 2 12 j ~ T ^ 720 a :2 “ 5040 ^ 


+ ... 


1 


Ao(a: + 1) = ( a: + - ) log a: - x + 


1 


1 


12a: 

1 


+ 


1 


360a:3 1260a:5 


( 9 ) 


1 1 

A 1 (a: + 1) = loga: H-^ + 

^ ^ 2a: 12a:2 120a:4 

. . X 111 1 

A-2(i + 1)--j + ^-^ + 


30a:^ 


I have included two ‘lower’ expansions which can be taken as part of the set. 
Explicitly, 

A-i + + 

A_2(a: + 1 )=p + ^ + -- - + ^-‘52 

where 5”^ and S 2 are the corresponding ‘Glaisher-Kinkelin-Bendersky (GKB)’ con¬ 
stants and have the (known) values. 


^i = -r'(i), 52 = C(2). 


We recognise in (9) some standard asymptotic (Stirling) series, e.g. for, 

Ao(a: + 1) = log , A_i(a: + 1) = ^ logr(a: + 1) = 'ip{x + 1) . 

and so on downwards to give the polygamma functions, up to a factor. 

As is visually clear, one can run up and down the right-hand sides in (9) by 
integration and differentiation. The exact relationship is given by Bendersky who 
derived it by brute force from (8). He found. 


^fc+i(^ -M) — (/c -f 1) 


Af.{x -t- 1) do: -|- 


j^(t>k+i {x + l)+xHf^ , 


( 10 ) 


where (pn{x) is a Bernoulli polynomial given in terms of the more usual polynomials, 
e.g. Norland, [5], B„(a:), by. 






which equals the sum of the nth powers of the hrst x — 1 integers. 
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Constants of integration can be considered to be absorbed into the GKB con¬ 
stants, Ljij., which are fonnd, in each case, by trial of varying x. 

It is important to note that, althongh derived nsing an asymptotic series, the 
recnrsion, (10), is valid generally. 

Another derivation of (10) is given later. 

The actnal n um erical valnes of the GKB constants, can be determined from 
the dehnition, (5) 

Lk = logr(«; + 1) - + 1) 

nsing (2) and a snitably chosen valne for the integer, w. This choice is linked to the 
necessary trnncation of the inhnite series in (8) and leads to an accnracy of at least 
29 places, e.g. for w = 100 and 20 terms of the snm retained. 

Bendersky also developes recnrsion relations for the which involve jnst con¬ 
vergent series. 

In order to hnd a more precise dehnition of log Tthan the asymptotic series, 
one starts from the observation that rQ(a: -|- 1) = r(a: -|- 1), in terms of the ordinary 
Gamma fnnction and takes this as the starting point of an npwards recnrsion based 
on (10). To take the simplest case, logr]^(a; -|- 1) has to involve the integral of 
log(r(a: -|- l)/\/^) and now, for exactness, a dehnite integral is nsed, 

Ytx + ll 

logr]^(a:-|-1) = / dx log-- \-X{x), 

Jo \/2 ti 

where X is to be fonnd. Becanse (10) has to hold np to a constant, C , 

X = (l)i{x -\-1) — C = ^x{x + 1) — C, 

which is determined by setting a; = 0, so that,^ 

C = X{tS) 

= 0 , 


and, hnally. 


logr]^(a:-I-1) = f dx logr{x + 1) + 4>i{x + 1) 

Jo 


X Lr 


— [ dx logr(a;-t-1)-I- 

Jo 


x{x -|- 1) 


X 


log 


( 11 ) 


The particular values logr^(l) = logr^(2) = 0 have been used as boundary conditions. 
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defining in terms of the standard F. (Setting x — 1 yields a known resnlt.) 

Barnes, [6] p.281, refers to this eqnation (or its eqnivalent) as Alexeiewsky’s 
theorem. 

Jnst to see how things ht together, I look at the next iteration according to 

( 10 ), 


r 11 

logF 2 (a: + 1) = 2 / dx logFi (x + 1) H— (p 2 {^ + 1) “I —^ ~ 2xLi 

Jo 2 6 

px px px 

= 2 / dx do;'log F(a:'+ 1) + 2 / dx[(l)i{x + 1) — xLq) + 

Jo Jo Jo 

+ + 1 ) + 

px px 2 

= 2 / dx dx'log F(a:'+ 1) H— (j) 2 {x + 1) — 2xL.^^ — x'^Lq . 

Jo Jo 2 

A continnation of this process prodnces Bendersky’s important general solntion. 


logFfc(x + l) = /c!4(^) + ^fc</'fc(^ + l)-V’fc(^), (12) 


with the dehnitions 


4(x) = / dx'4_i(x'), 

Jo 



h{x) = logr(a:+ 1) 


(13) 


Bendersky takes the convergent (12) as the dehnition of the generalised Gamma 
fnnctions for all x and any integer k and nses it systematically to obtain their basic 
properties snch as (7), enconntered earlier. It wonld be ont of place to give a 
snmmary here. 


4. Earlier history 

The qnantities logFj!.(x) in (2) were hrst considered by Kinkelin who concen¬ 
trated on jnst F 2 (x) in [7],^ and derived its properties ab initio from the expression 
in terms of the ordinary Gamma fnnction, 

logF4x) = y' logF(t) dt + - ^x log27r (14) 

^ Kinkelin’s work is dated 1856, but was published in 1860. I have not been able to access his 
earlier paper. 
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which he obtained using Raabe’s formula and the fact that it reduces to the sum¬ 
mation, (2), for X integral. 

He computes what is essentially the constant not by the trial method but 
from derived convergent series. These series can also be found in Bendersky. 

Glaisher, [8], gave a numerical treatment of using the trial method. In 
this, however, he was forestalled by the earlier work of Jeffery ^ in a known paper 
[10], written in 1860. In a second, neglected paper, [3], written in 1862, Jeffery 
gives what seems to be the earliest calculation of the second GKB constant, L 2 - 
Indeed he describes, inter alia, the construction of the complete general recursion 
system, (10), and I have taken equation (9) from his paper. It is clear he could have 
computed as many of the Lf^ as desired and so I propose that GKB be extended to 
GKBJ. 

It is worth noting, with Jeffery, that the coefficient of log a: in -|- 1) equals 
or 5i-|-fc(ir)/(l -|- k), although he does not seem to prove this in generality. 
Jeffery’s approach is worth re-exposure. It does not overtly use the Euler-Maclaurin 
formula. 

5. Jeffery’s treatment 

For shortness, I dehne 


Mx) = logrfc(a: + 1) 

From the original expression (2) algebra gives, 

= ix + 1)'" log{x + 1) , 

which dehnes a system of equations. 

Extending x into the reals gives 

DAvi^{x) = ADvj^{x) = k{x + 1)^“^ log(a: -|- 1) -|- (x -|- 1)^“^ 
= kAv^_i{x) + [x + 1)^“^, 

so that, by summation, 

Dv^{x) = kA~^Av^_]^{x) -|- 1)^“^ -|- W]^{k) 

= kv^_^{x) + + 1) + w{k) 

^ A very brief biography can be found in [9]. 


(15) 


( 16 ) 
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The final periodic constant w{k) is given by setting x = 0 whnch yields 

zu{k)=DvkiO), /c = 2,3,... 

• 07 ( 1 ) = D Vi{0) — 1 


( 17 ) 


since ^^^(0) = 0 = I will nse (17) to find n;.(0) rather than the w{k). 

Jeffery, [10], gives only the nnmerical valne iu{2) — —0.2475089541... which 
is 1/4 — 2Li, in terms of the GKBJ constant for k = 1. This can be shown in the 
following somewhat lengthy way. How Jeffery works ont the valne is not clear to 
me. 

Integrating (16), 

fX 

Vj.{x) — C{k) = k / dx Vj^_i{x) + {x + 1) — Bf.x) + zu{k) X (18) 
Jo ^ 

where C{k) = 0, again on setting x = 0. 

As a trivial check, set k = 1. Then 

fX 

Vi{x) = / dxvQ^x) + {4>i{x + 1) — B^x^ + w{l) X 

Jo 

1 7 7 +1) 7 717 I 7 

= dxvo{x) + - - — + {w{l) +-) X . 

Setting X to zero prodnces nothing while x — I gives, nsing Raabe’s formnla. 


0 = / dxvQ^x)+ w{l) 
.In 2, 


= log + ti7(l) + ^ 


therefore the relation to the GKBJ constant, Lq, is 


^(l) = -^0 - i = - log ^ - I ■ 


Then 


Vi{x) = / dxvQ{x){(l)i{x+ 1) — B^x)-\-w{l) X 

Jo 

1 7 7 + 1) 1 /TT- 

= / dxvQ[x)^ ---iogv 27 ra:. 


(19) 


which is the same as (11), or (14). 

Incidentally, from (17) , one sees for that the series for Vi in ascending powers 
of X begins with the term (l/2 — log \/2jr) x. 



Relation (18) can be iterated. Thns 


Vkix)=kJ^ dx' dx"v,^_ 2 {x”) + -^^{(j)k-iix' + l)-B^_^x') + 

+ w{k — 1) x' + — {4>k{x + 1) — B^x) + w{k) X 

fX r-x' ^ 

= k{k-l)J dx'J dx"v^_ 2 ix") +-^^{(j),^{x + l)-B,^x) 

- ^ ~ ^ 1 + ^)-Bkx)+ w{k) X . 

(20 

The simplest, non-trivial, case is /c = 2, when, 


oX px 


V2{x)=2 / dx' / dx" Vq{x") + {(j) 2 {x + 1) — B 2 X^ 

Jo Jo 

— BiX^ + tu{l) H— {(f>2ix + 1) — B 2 x) + w{2) X . 

2 

Setting X — 1 enables an expression for the summation constant, w(2)^ to be found, 


0 = 2 / dx' / dx" Vq{x") + ((/>2(2) — B 2 ) 

Jo Jo 

- B^+ w{l) + ^((/)2(2) - B 2 ) + w{2) 

= 2 / dx{l — x) logr(a: + 1) + { 4 > 2 i‘^) ~ -^ 2 ) 

Jo 

- B^+ zu(l) + ^ ((f> 2 i 2 ) - B 2 ) + tu{2) 

7 

= 2 dx {1 — x) log r(a: + 1) + - + '07(1) + - 07 ( 2 ). 

Jo 4 

Rewrite this in terms of r(a:), 

3 7 

= —2 dxx log r(a:) — - + 2 log + 7 + ^(1) + '07(2) 
Jo 2 4 

/•i 3 

=—2 dx a; log r(a;) — - — • 07 ( 1 ) +- 07 ( 2 ). 

7n 4 


In order now to relate the constant of integration, t<7(2), to the GKBJ asymptotic 
constant, L^, I nse Kinkelin’s general approach. 
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The derivation is based on the classic expression, or dehnition, of the gamma 
fnnction as an inhnite prodnct. This gives, 

M —1 

logr(a:) = logw! + (x — 1) logrc — log(a: + t), w ^ oo . (23) 

i=Q 

Also needed is the asymptotic limit, assnmed known, of the factorial (Stirling), 


log w\ ~ —w + w log w + - log w + Lq . 
The relevant integral in (22) is 


„l ^ U)-l „l 

2 / dxa; logr(a:) = logrc!-log re —2 / dx x \og{x-\-i), (25) 

Jo 3 Jo 


and so one needs. 


W — l„l W — lx 

2 ^ / dxx log(a: + i) = ^ (log i — log(t + 1)) + log(t + 1) + t — 

,_n Jo ,_n V 


Rorganising, 


( i‘^[logi — log(t + 1)) + log(t + 1) j = ^ logt — ^(*^ — 2t + 1 — 1) logt 
i=o ^ J i=l i=l 

W 

= 2i log i — w‘^ log w , 


giving, for this piece. 


w-l „1 


2 / dx X \og{x + i) = ''^‘^2i\ogi — w'^ \ogw-\ - 

,-n Jo 2 


Combining this expression with (25) and (24) I hnd. 


2 / dxx logr(a:) = 


w log w + 


logre 1 


6 2 


- --07(1) - ^2tlogt + M;^logM; - — . 
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Then (22) becomes 


3 

4 


tu{l) + w{2) 


w log w + 


log w 
6 



W 

(1) - ^2ilogi + 

i=l 


w‘^ log W — 


W 


(27) 


or 


^ -1 f ^ ^ / 2 ^^ 1 

2_^ilogi = I — + -(w + m; ) 1 logM;+ - - ^ - 

i=l ^ ' 


-w(2) 


showing that the constant, is, by dehnition. 




The initial conditions (17) also then show that the power series expansion of 
V 2 {x) = logr 2 (a: + 1) starts with (1/4 — 2L-^)x as stated, nnmerically, by Jeffery, 
[3] §43. 


Kinkelin derives a constant, to, as a convergent series, from mnltiplication prop¬ 
erties of and then shows that this is the constant in the asymptotic form. It can 
be fonnd from the information already obtained if I assnme Kinkelin’s eqnn.(22“) 
which is. 


log TO 


1 

~6 

2 

~3 

1 

12 


+ log\/^—2 / dxx\ogT{x) 


-^( 1 ) + ^ +^( 1 ) - 

-to( 2 ) = 2 Li - ^ ~ 

D 


- to (2) 

0.33084228740 


and agrees with Kinkelin’s valne, obtained by series snmmation. 

Snch a case by case approach, directly applied, is not algebraically efficient for 
the higher iterations. As two hnal examples I hnd that. 


to(3) = - 3 L 2 ~ -0.091345371176, to(4) 


- - 4Lo ~ 0.013180972097. 

72 ^ 
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6. The general solution. Some integrals 

I briefly return to Bendersky’s general solution, (12), which he takes as the 
deflnition of the function. The ingredients are in (13), the main part of which is the 
nested integral which can be reduced to a single (fractional) integral of log T, 

1 

4(^) = 1); dt{x-t)’^-^ logT{t + l). (28) 

By giving x particular values, expressions for this integral can be found. For 
example if a; = 1 the left-hand side of (12) vanishes and so, 

k f + = 

Jo 



(29) 


As a check take k = 2. Then simple arithmetic yields, 

2 dt (1 - t) log m = -C'(0) - 2C'(-1) + ^, 

which agrees with the value given in [11] p.674. In fact equn.(109) of [11] contains 
(29). Other integration ranges can be accommodated. 

I note that Ij^, (28), is directly equivalent to the polygamma function of negative 
order, discussed by, e.g., Adamchik, [4]. See also Espinosa and Moll, [12]. 

More precisely, just changing r(t + 1) to r(t), 

4(4 dt {x - t)^~^ logr(t) + ^ dt [x - \ogt 

—^ [ ((x — — x'^)t~^dt + x'^ logx 

k'- Jo 

= (X) + ^ + log:,]) , 

Furthermore, Proposition 2 in [4] for the polygamma function is equivalent to 
the general solution (12) together with the basic relation, (1). 
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7. Comments 


Unless I am missing something, Jeffery’s method of section 5 does not easily 
give the general form for Vf^. Furthermore, in the details of the algebra, cancellations 
occur which make w{k) a function of Lj^_i only, which ought to be derivable directly. 

By contrast, Bendersky’s approach is a top down one in that a general recnrsion 
is found, somewhat out of the blue, for Af. = Vj^ — Lj^. This leads, essentially by 
induction, to the form of Vj^ for all k, (12), which is then taken as the definition of 
Vf,{x) for any x. Equation (12) straightaway shows that, 

— kLj^_i 

as confirmed above in particular cases. A derivation of this expression that does 
not depend on the explicit form of the general solution for Vj^ would be welcomed. 

8. Alternative approach 

I mentioned that looking at the basic equation (1) as the definition of rj!.(a:) 
allows a rapid derivation of its properties and relations. As an example, I now derive 
the recursion (10) which follows almost immediately from the standard formula, 
used in [11] for a related purpose (see also [13] for a generalisation), 

^C(s,«') =-sC(s+ !,«'), (30) 

so that, differentiating with respect to s, 

^C'(s, w) = -C{s + l,w) - s C{s + l,w). 

Then, setting s = —k with k a positive integer, and integrating, gives the result, 

(^ + 1) ^ (c'(-^, t) - dt = C'{-k - 1, x) - C'{-k - 1). 

where I have used ('{-k, 0) = ('(-k), [11], Appx.C. 

This recursion is equivalent to Bendersky’s, (10), as a few lines of algebra, util¬ 
ising (3), confirms. Kurokawa and Oshiai, [14], unaware of the work of Bendersky, 
use this approach and also generalise it to Barnes C^functions. I note that their 
definition of the generalised gamma function (‘higher depth’ gamma function) does 
not include the (constant) second term on the left-hand side of (1). 
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As a minor point, one of the integrations can be done using the well known 
relation, 


({—k — l,x) — ({—k — l) = {k + l) f dt({—k,t), k>0, 

Jo 


which follows trivially from (30). It expresses a standard relation between Bernoulli 
polynomials. 

Continuing this theme of rewriting, Adamchik’s form of the GKBJ constants 
is, 

L, = -C{-k)-H,a-k). (31) 


But I take this no further since, I will give, at another time, an extended analysis 
of this approach and its generalisations as in [11], [15] and elsewhere. 


9. Summary 

I have shown that Glaisher’s numerics and Bendersky’s recursive approach 
to the generalised gamma functions were both anticipated by Jeffery whose work 
appeared a couple of years after Kinkelin’s paper, of which Jeffery seems unaware. 

I have advocated the computation of the derivative of the Hurwitz (^-function 
at the negative integers via the generalised gamma function, rather then vice versa. 
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